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INTRODUCTION 

>  Problems  of  elastic  wave  scattering  by  surface-breaking  and  near-surface  cracks 
are  of  considerable  current  interest  for  ultrasonic  nondestructive  evaluation.  Ultra¬ 
sonic  scattering  by  .planar  cracks  near  or  at  the  free  surface  of  a  semi-infinite  elastic 
homogeneous  medium\has  been  studied  theoretically  by  many  authors.  References 
to  recent  papers  on  this  subejct  can  be  found  in  [1],  [2]  and  [3].  Some  experimental 
works  on  surface-breaking  normal  planar  cracks  have  also  appeared  ([4]  -  [6]). 

Ultrasonic  scattering  by^urface-breaking  planar  and  branched  cracks  of  arbitrary 
orientation  is  the  subject  of  this  investigation.  YTo  our  knowledge  this  problem  has  not 
received  much  attention  in  the  literature.  An  approximate  solution  that  is  valid  at  low 
frequencies  was  presented  in  [7]  for  SH  wave  diffraction  by  a  canted  surface-breaking 
planar  crack.  Subsequently,  a  hybrid  finite  element  and  eigenfunction  technique  was 
used  in  [8]  to  study  SH  wave  diffraction  by  planar  surface-breaking  canted  crack. 

— ■' '  /  In  this  paper  we  use^the-sarae  hybrid  t  .chnique  as  in  [8-fHo  study  the  scattering 
of  in-plane  body  and  surface  waves  by  canted  planar  and  normal  surface-breaking 
branched  cracks.  We  focus  our  attention  to  the  near-field.  Numerical  results  are 
presented  for  the  vertical  surface  displacement  amplitudes  near  the  base  of  the  crack.  — ^  j 

FORMULATION  AND  SOLUTION 

Consider  a  homogeneous,  isotropic,  and  linearly  elastic  medium  with  a  near¬ 
surface  inhomogeneity  of  arbitrary  properties  and  shape  as  shown  in  Figure  1.  Assume 
the  displacement  u(z,y,t)  at  a  point  P  to  be  time-harmonic  of  the  form  u(x,y)e“‘wt, 


Fig.  1  Geometry  of  a  near-surface  inhomogeneity 

where  u  is  the  circular  frequency.  Then  u  satisfies  the  equation  of  motion  in  y  >  0 
(at  points  not  on  the  crack) 

/iV2u  +  (A  +  ft)  VV  •  u  =  —  pw2u  (1) 

where  \,fi  are  Lamd  constants,  p  the  mass  density  and  the  factor  e~iut  has  been 
dropped. 

The  solution  of  (1)  can  be  expressed  in  terms  of  the  longitudinal  and  shear  wave 
potentials,  <f>  and  in  the  form 

tt  =  V^  +  Vx  (ipe,)  (2) 

Furthermore,  in  a  homogeneous  half-space,  <P  and  rp  can  be  expressed  in  infinite  series 
of  multipolar  potentials  as  [9], 

OO 

^  =  ^2  [an<Pn  + 

»s-oo 

i>=  J2  (®»^  +  Wi)  (3) 

ns— co 

where  expressions  for  and  ipX  can  be  found  in  {9).  The  coefficients  an,  bn 

are  found  by  satisfying  the  appropriate  boundary  conditions. 

The  representation  (3)  is  not  useful  for  satisfying  the  boundary  conditions  on  the 
crack  surface.  For  this  reason,  a  different  representation  is  needed  in  this  near-field 
region.  In  this  paper,  the  region  inside  the  fictitious  boundary  B  (Fig.  1)  is  divided 
into  finite  elements  having  N[  number  of  interior  nodes  and  Np  number  of  boundary 
nodes. 

For  the  finite  element  representation  in  region  II,  the  energy  functional  is  taken 

to  be 

F  =  \  [  f  (a  •  «*  -  pu2u  •  u*)dxdy  -  ^  f  {Pd-ud*  +  Pd*  •  uB)ds  (4) 
2  J  Jr,,  2  JB 


where  “  *  *  indicates  complex  conjugate  and  a  and  e  are  column  vectors  defined  as 


T 

O  —  {0}  =  {p tx  >  &yyi  &xy  ) 

(5) 

t  —  {*}  =  (e**»  <yy»  t-xy) 

(6) 

Superscript  uTn  denotes  transpose.  The  Pq  and  UB  represent  the  traction  and  dis¬ 
placement  on  B,  respectively. 

It  is  assumed  that  the  displacement  field  within  the  jth  element  is  represented  in 
terms  of  the  shape  functions  L}(x,y)  and  elemental  nodal  displacements  {<7*}  as 

£  L>1"  (7) 

y=i 

where  each  q‘  has  two  components  uX]  and  uyl-  along  the  x  and  y  directions,  respec¬ 
tively.  The  N,  represents  the  number  of  nodes  in  each  element. 

The  and  are  computed  by  substituting  Eq.  (7)  into  strain-displacement 
relations  and  these,  in  turn,  into  the  stress-strain  relations.  Using  these  in  Eq.  (4), 
we  get 

P  —  q]T  S I  iQ  I  +  9/T,S’/b?B  +  9bT  ^b/9/ 

+?bT  ^bb?b  “  Q*b  “  P b'^9b  (8) 

in  which  qt  =  q^2\qB  =  qg^KPg^1^  -  Pg^  the  elemental  impedance  matrices 
Sij  are  defined  as 

l«1  =  /  -  ,*,\l\T\L\)dziy.  (9) 

In  Eq.  (9), 


Note  that  [L]  is  a  2  x  2 N,  matrix. 

For  an  isotropic  material  [D]  is  given  by 

A«  +  2  Ht  A*  0 

[D\  =  A.  A.  +  2m.  0 

0  0  M*. 

where  A,  and  n,  are  the  Lamp’s  constant. 

To  find  the  constants  an,bn  appearing  in  (3)  and  the  nodal  displacements  in 
region  II,  it  is  necessary  to  use  the  continuity  of  displacement  and  traction  on  B. 
This  is  discussed  in  the  following. 

The  incident  displacement  fields  will  be  assumed  to  arise  from  the  incident  plane 
P  and  SV  waves,  and  their  reflections  from  the  free  surface  y  =  0.  The  case  of 
incident  Rayleigh  waves  will  also  be  considered. 


Let  us  suppose  that  in  the  absence  of  the  crack  the  free  neid  is  the  sum  of  the 
incident  and  reflected  fields,  that  is 

u?> = «<*> + «;.r|  (j  =  i,2).  no) 

For  the  Rayleigh  wave  u*01  is  the  associated  displacement. 

The  total  field  outside  B  then  is 

uy  =  uf>  +  u'0)  0  =  1.2).  (11) 


where  uj5*  is  given  by  (2)  and  (3). 

Using  (2)  and  (3),  the  displacements  at  the  nodes  on  B  can  be  written  as 

{«if}  =  (C]{o}  (12) 

where  [G]  is  a  2 Nb  x  2 Nb  matrix  formulated  in  [10]  and  vector  {a}  is 

T 

=  [®  1.  •••>  &Nb  >  ^1 » •••>  bftg  ] 

Similarly,  using  (2)  and  (3)  in  the  stress  strain  relation,  the  traction  at  the  nodes 
on  B  can  be  expressed  in  the  form 

=  mw  (i3) 

where  [/*]  is  also  a  2 Ng  x  2 Ng  matrix  defined  in  [10]. 

To  express  {o^}  in  terms  of  we  use  expression  for  the  virtual  work 

done  on  the  boundary  B,  which  is 


Sir  =  f  {Vl1)}aW1)}«*r  (14) 

Jb 

where  superscript  (1)  denotes  the  total  field  in  region  I  (outside  B). 

Because  of  the  continuity  of  displacements  and  traction  on  B ,  we  have 


(15) 

(16) 

where  superscript  (2)  denotes  the  total  field  in  region  II. 

Substituting  (12),  (13),  (15),  and  (16)  in  Eq.  (14),  and  noting  that  6q ^  * 
we  obtain  from  eq.  (14) 

II 

o> 

5  *3 

where  is  given  by 

(rt11)  =  !*!{«}  +  (<') 

(18) 

Here 


and 


1*1  =  /  [G'iT[F)<ir 

Jb 

{*>£”}  =  /  (c,|T{40)>Jr. 

Jb 


(19) 

(20) 


Substituting  Eq.  (12)  into  Eq.  (8)  and  taking  the  variation,  we  obtain  a  set  of 
simultaneous  equations  which  may  be  written  in  matrix  form  as 


Sji  SibGi  qi  _  —StbQb^ 

G-tSJ„  g-ts°b  J  [  J  L  -G’tSb „$>  +  P‘Ba>  ■ 


(21) 


Using  Eqs.  (18)  and  (21),  we  obtain 


[G*T(Sfl* -SJbS^SjbG*  -  £]{o}  = 

-G-t{SBB  ~  SjBSrIlSIB)g^)  +  Pj,0*  (22) 

In  Eq.  (22),  the  generalized  coordinates  {a}  are  the  only  unknowns.  Therefore,  {a} 
can  be  evaluated.  Once  {a}  are  known,  the  near  and  far  displacement  and  stress  fields 
can  be  determined. 

NUMERICAL  RESULTS  AND  DISCUSSION 

In  this  paper,  the  boundary  B  enclosing  the  interior  region  is  not  a  complete 
circle,  and  so  the  potentials  and  cannot  be  expanded  in  cylindrical 

wave  functions  as  was  done  in  [lj.  So  the  integrals  giving  these  potentials  and  their 
derivatives  were  evaluated  numerically  for  every  node  on  B.  The  details  are  discussed 
in  [11]. 

The  hybrid  method  is  employed  to  study  scattering  by  P,  SV  and  Rayleigh  waves 
by  three  types  of  surface  breaking  cracks:  a  vertical  crack  (Fig.  2  with  a  =  90°),  a 
45°  inclined  crack  (Fig.  2  with  a  =  45°),  and  a  vertical  branched  (Y)  crack  (Fig.  3). 


Fig.  2.  Geometry  of  a  planar  crack 


Fig.  3.  A  branched  crack 


Stress  intensity  factors  at  the  tips  of  the  cracks  were  calculated  and  for  the 
particular  case  of  a  planar  surface-breaking  normal  crack  they  were  found  to  agree 
well  with  the  results  of  [12]. 

The  surface  displacements  at  y  =  0  are  calculated  by  using  (3)  in  (2)  after  {a} 
are  calculated.  Normalized  values  of  UyS^  are  presented  in  Figures  4-8.  For  each  type 
of  crack  mentioned  above  five  cases  of  incident  waves  were  considered:  plane  P  wave 
incident  at  0°  and  45°,  plane  SV  wave  incident  at  0°  and  45®,  and  finally  Rayleigh 
wave.  Some  representative  results  are  shown  here. 

Figures  4  and  5  show  the  scattered  vertical  surface  displacement  amplitudes  for 
a  Rayleigh  wave  incident  from  the  left  on  a  normal  planar  and  branched  crack.  It  is 
seen  that  there  are  large  differences  in  the  forward  direction  between  the  two  cases 
as  the  frequency  becomes  large.  In  the  backward  direction,  however,  the  differences 
are  not  very  significant.  Figures  6  and  7  show  the  results  for  an  incident  SV  wave 
moving  vertically  as  well  as  at  45°  to  the  vertical.  Large  differences  are  found  for 
vertical  incidence,  but  not  in  the  other  case.  Finally,  in  Figure  8  is  shown  the  case  of 
a  Rayleigh  wave  incident  on  a  canted  crack.  This  figure  is  to  be  contrasted  with  Fig. 
4.  The  large  contrast  shown  clearly  distinguishes  a  canted  crack  from  a  normal  crack. 

CONCLUSION 

^  Model  calculations  of  elastic  wave  scattering  by  surface-breaking  planar  and  non- 
planar  cracks  have  been  presented.  These  calculation^  show  that  near-field  surface 
displacements  due  to  scattering  by  planar  and  branched  cracks  are  quite  different  even 
when  the  branches  are  small.  Also,  it  is  found  that  signatures  of  normal  and  canted 
cracks  are  very  dissimilar.  These  characteristic  differences  can  be  used  to  discriminate 
between  the  various  cases. 

% 

Although  the  results  presented,  here  are  for  homogeneous  medium,  the  technique 
can  be  generalized  to  study  cracks  in  a  composite  medium.  These  are  presently  under 
investigation  and  will  be  reported  later.  -/ 
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Fig.  4.  Scattered  vertical  surface  displacement  amplitude  due  to  normal  planar  crack 
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Fig.  5.  Scattered  vertical  surface  displacement  amplitude  due  to  normal  branched  crack 
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Fig.  8.  Scattered  vertical  surface  displacement  due  to  45°  inclined  crack 
for  incident  R  wave 


